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ABSTRACT 


The Naval Postgraduate School has an ongoing effort to study 
geomagnetic noise and micropulsations in the ULF frequency range 
(.05 4 £<10 Hz). Data is collected by three orthogonally 
mounted coils at a remote land site and telemetered to the 
laboratory for computer analysis. 

To isolate data containing a micropulsation event, time 
series plots of the magnetic field were generated. The 
development of a double running average routine made possible 
the isolation of micropulsations in large data sets. A type 
PC-3 micropulsation was found and the coherence, ellipticity and 
polarization properties were determined as follows: Coherence = 
0.99, Degree of Polarizaticn = 9.99 and the micropulsation was 
elliptically polarized. 

Power spectral density (PSD) plots summarizing data segments 
about two hours long occasionally contained structures found to 
be artificial and not representative of natural phenomena in the 
geomagnetic field. Methods to avoid this anomalous behavior in 


PSD plots are suggested. 
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This thesis is part of an ongoing effort at the Naval 
Postgraduate School to analyze ULF geomagnetic noise and 
micropulsations. The results of the Naval Postgraduate School 
studies could impact communications systems or systems in which 
geomagnetic noise and/or micropulsations introduce operational 
ett tculties. 

An objective of this thesis was to develop the software 
necessary for obtaining computer generated plots of geomagnetic 
field versus time. Then, with the plots, sections of data 
containing obvious geomagnetic micropulsations were located. 
Those sections were analyzed to determine the power spectral 
density (PSD), coherence, degree of polarization and 
ellipticity of the micropulsation using already developed 
software [Ref. 1]. 

After analyzing some of the plots (PSD, coherence, etc.) 
oroduced by the software documented in [Ref. 1], it was 
demonstrated that problems existed in the software. New 
software was developed to obtain times series raw voltage and 
geomagnetic field plots. Validation of the software required 
an on-site experiment. Documented within this thesis is the 
software to produce time series voltage and geomagnetic field 
plots and the questions the validation experiment raised 


concerning the sensing system and software. 





Using the newly developed software, a micropulsation event 
was located. A double running average software routine 
accentuated the micropulsation to permit the determination of 
its frequency and amplitude. The micropulsation event was 


further analyzed to obtain information on coherence, degree of 


PolariZation and ellipticity. 





If. BACKGROUND 


Pee COLLECTION SYSTEM 

The Naval Postgraduate School has two operating geomagnetic 
sensing systems. The land site is located at the La Mesa 
Village Housing Area and has three orthogonally mounted 
induction coils. The second operating system contains two 
orthogonally mounted induction coils that can be placed on the 
bottom of Monterey Bay. When both systems are operating, the 
data is collected simultaneously. The ground work has been laid 
for a third system located at Chew's Ridge. Three orthogonal 
ficmelinmlcaneous measurements can be taken only between two 
Sites since the School has only five induction coils and a sixth 
just calibrated. This thesis used data from the La Mesa Village 
Site and the software associated with it. 


An induction coil senses the geomagnetic field fluctuations 


ab 
ct 


deseribpes the relation between the induced emf and the time 


based on Faraday's Law of induction, induced emf =-N which 
maieewor Change of magnetic flux through the coil. N is the 
Mmumeer Of turns in the coil. 

The voltage induced in the coils is sampled 32 or 64 
times/second and is amplified approximately one million times. 
The voltage data undergoes pulse-code-modulation (PCM) for 
noiseless VHF radio link transmission from the La Mesa 


Village site, to the recording system in the Geophysics Signal 


10 





Pe@ecessing Laboratory located in Spanagel 531. Data reduction 
on the Naval Postgraduate School's IBM 3033 computer follows 
digitization of the PCM data. A block diagram of the ee is 
provided in Figure 2.1. A more detailed description of the 
system 1s given in Reference 2. 

1. System Calibration 

For purposes of documentation, the calibration of the 
third coil for either the Chew's Ridge or Monterey Bay system 
is being included in this thesis. The general method of 
induction coil calibration will be covered here. For a more 
indepth explanation, see Reference l. 

The sensing coil and its associated electronics are 
calibrated using a Helmholtz coil apparatus to establish a 
uniform magnetic field. By placing the sensing coil into the 
field produced by the Helmholtz coil, one can establish a 
relationship between the Helmholtz coil magnetic field magnitude 
and frequency and the coil system voltage. Once the response of 
the coil system is known for the frequency band of interest, it 
can be included in the system transfer function used in data 
reduction. 

The experiment to determine the sixth coil system 
transfer function was done at the La Mesa Village site. A 
Wavetek Model 142 signal generator supplied a sinusoidal current 
tems 22 meter diameter, 0.61 meter high Helmholtz coil. The 
current was measured across a 996 ohm resistor in series with 


the Helmholtz coil. A Hewlett Packard HP-3582A spectrum 
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analyzer measured the output voltage of the coil system as a 
BeMeGelon Of frequency. 

The experiment measured the response of the ee system 
in the frequency band .05 hertz to 20 hertz for applied fields 
of 0.02, 0.2 and 1.0 nanoteslas. The coil system response is 
Shown in Figure 2.2. The transfer function of this coil system 
corresponds closely to those of coil systems previously 
calibrated. The only noticable difference is that the 0.02 
nanotesla response at 15 hertz is approximately 15 volts/ 
nanotesla higher than those measured before [Ref. 1]. It 
should be noted that this difference will not influence data 
analysis in the frequency range of interest. The transfer 


function algorithm is listed in Table 2.1. 


pe. PREVIOUS SOFTWARE 

Software has been develoned (see References 1 and 3), to 
extract power spectral densities for each coil, coherence 
between the coils, degree of polarization in the three measure- 
ment planes, ellipticity in the three measurement planes and a 
variety of Stokes parameters. The following is a basic 
description of the computer code that produces the above 
mentioned plots. First, the voltage data is read off the 
digital tape using a subroutine called RD provided by Dr. Tim 
Stanton of the Oceanography Department. A parameter ISEC, 
representing the number of seconds one wishes to advance the 


tape, is frequently used. ISEC establishes the number of 


13 





Sf 


uoT}JOUNT AAOySsuerzy, TTOD 


a YIH) ASNINGAYS 


nn nm es ET) ee on capi 


LUCO* 


LUT 


LU ° 


Wad Saher esl 


a I 


NOILANNA SAASNHYI 


CINSA? 


14 





system Transfer Function Algorithm 


Table 2.1. 
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times RD is called and the data not stored. Previously, tape : 
advances of approximately 30 seconds or ISEC = 30 were -used. 
Once the tape is advanced the desired amount, data analysis can 
begin. Each digital tape contains about 90 minutes or 1.728 x 
10° pieces of data per axis resulting from the sampling rate of 
32 samples/second. Because of the large amount of data, the 
analysis has to be accomplished in blocks. Data analysis takes 
place on 256 seconds or 8192 frames of data at a time. RD is 
called 8192 times and places the voltages for each coil in an 
array. The voltages are integer values from 0 to 4096 
representing -5 volts to +5 volts. Ideally the value 2048 
represents O volts. The integers between O-4096 are normalized 
to + 5 volts and placed in new complex arrays. These arrays 
containing 8192 pieces of data are Fourier transformed to the 
frequency domain where the system transfer function is applied. 
After applying the transfer function, the calculations begin 
Bamerecetermination Of PSD's, Coherence, ellipticity, etc. for 
the first block. After they are complete, the next block of 
8192 pieces of data starts through the program. Once the 
desired amount of data has been processed, a NONIMSL subroutine 
called DRAWP plots the data. For a more complete description 


of this software, see Reference 1, Reference 3 and 


Reference 4. 
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Se aNAGbYSIS OF INITIAL DATA PLOTS 
Many of the PSD's, coherences, ellipticities and degree of 


polarization plots produced by the software of Reference 1 and 


Reference 3 have unexpected characteristics. The PSD's show 
"humps" that are shown in Figures 2.3-2.5. These "humps" are 
ieememasacteristic In “normal PSDs. The coherence plots show 


a coherence of essentially one degrading to hash in the 
Frequencies greater than 2 hertz, Figures 2.6-2.8. We see a 
Similar behavior in the plots of degree of polarization. These 
plots indicate a very highly polarized field from 0.02 hertz to 
memimersez. Figures 2.9=-2.11. The elliptici#y plots show a very 
linearly polarized field for these frequencies in Figures 
Ze12—-2.14. 

After examining the time series magnetic field plots, we 
noted a section of very large magnitude fluctuations (about one 
to three orders of magnitude greater than the general back- 
ground) at the beginning of the digital tapes Figures 2.15-2.17. 
We found that if the digital tape was advanced past these large 
fluctuations the output (PSD, coherence, etc.) appeared 
“normal", Figures 2.18-2.29. If the program was executed on 
just the section containing the large fluctuations, the PSD's 
had the "humped" characteristic, the coherences were one, the 
degree of polarization was very high and the ellipticity was 
linear, Figures 2.30-2.41. This demonstrates that these large 
fluctations dominate the output of an analysis in the frequency 


domain. A tape advance of 300 seconds or ISEC = 300 is 
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recommended at the beginning of the programs producing tne above 
mentioned plots. Time series plots should also be made to 
provide a further check. 

An explanation of the source of the large fluctuations is 
that the PCM to digital conversion begins before the actual PCM 
data starts on the tape. The decoding equipment considers the 
noise on the section of tape before the PCM data starts as 
Signal. 

The above analysis clearly proves that the "hump" like 
structure in the PSD's and the "nice" behavior observed in the 
corresponding coherence, ellipticity and degree of polarization 
are artifical and do not represent any real phenomena in the 
Geemagnetic field. Thus, the data in Figures 2.3-2.14 and 
Figures 2.30-2.41 does not appear to be reliable. Unfortu- 
Paieeny, this artificially induced, and therefore unphysical, 
behavior of the polarization has been published recently 
[Ref. 7]. The conclusion found in Reference 7 that 

Geomagnetic fluctuations in land are well polarized 

below 1.0 Hz and that the polarization is quite 
frequently, though not consistently, linear. The 
pisclatlonm under the sea is similar. Furthermore, 
a land-sea coherence study in the horizontal plane 
showed very high coherence below 1.0 Hz. 

are not substantiated by the physical data. Plots of the 


Paeeereorahiy induced behavior are bublished in References 1, 2, 


ama 3. 
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The 4 HZ peak present in the PSD's shown in Figuiiecee a= 
oe0 1S the result of 60 Hz aliasing. 60 Hz magnetic Peds 
generated by local power lines exist at the measurement site. 
Digitization at a 32 Hz or 64 Hz rate of the analog voltages 
from the sensing coils shifts the true frequency of the power 
beee signals down at 4 Hz. The 60 Hz aliasing’ results in the 
anomalis behavior at 4 Hz seen in the coherence, degree of 
polarization and ellipticity plots shown in Figures 2.21-2.29. 

Another peak is present at about 1.2 Hz for the PSD's in 
Poawecs 2.18-2.20. For this data set, the 1.2 Hz signal 
Semrelates with a 1l.2 HZ signal identified in the time series 
Soe. (ne Origin of the 1.2 Hz signal is unknown. The 4 Hz and 
ies peaks are of little concern in this study since both are 
well above the frequency of the micropulsation analyzed. 

Another "suspicious" feature of the PSD's is the large 
amount of hash appearing from 2 to 10 hertz. We will gain 
further insight into this problem in analyzing the time series 
voltage data. 

From the analysis of the abnormal behavior in the PSD's, 
coherence, etc., future coil data sets should receive a 
balanced treatment between the time domain and the frequency 


domain to avoid misinterpretation. 
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Iii. TIME SERIES DATA 


A. TIME SERIES VOLTAGE DATA 
1. Voltage Software 

The importance of time series data as a check became 
evident when we uncovered the problem with the initial PSD's, 
coherences, etc. Computer code was developed to produce both 
unsmoothed and smoothed times series voltage plots. The reason 
for producing both smoothed and unsmoothed voltage plots will 
become apparent in the time series voltage analysis section. 
Time series voltage plots could provide a means of monitoring 
the data for system "glitches" which might later produce 
erroneous PSD's, etc., and suspicous features on time series 
magnetic field plots. 

The Fortran computer orogram that produces the voltage 
plots is compiled and run on the Naval Postgraduate School IBM 
3033 VM or "batch" system. Both the unsmoothed voltage program 
(VOLTR) and the smoothed voltage program (VOLTS) require 2.048 
megabytes of core memory and approximately 15 minutes of 
Semecale processor unit (CPU) time to run. These figures 
provide a comfortable margin to insure that the program does 
not "bomb" due to lack of core or run time. 

The design of both VOLTR and VOLTS is the same with the 
exception that a smoothing algorithm is applied to the data in 
the case of VOLTS. Both programs have incorporated in them the 


digital tape advance package using ISEC as the number of 


a0 





seconds one wishes to advance the Gape.  Pataweiswreac Ofrt the 
digital tape in blocks of 8192 frames by the subroutine RD 
referenced earlier. As before, this data is in integer form 
and represents voltages between +5 volts by numbers between 0 
and 4096. The integer representing 0 volts is ideally 2048, 
however, I found that the "real" zero values wander. The 
reason for this "wander" is not understood. For good results, 


the following values worked well for the respective coils. 


X-coil : zero = 1966 
Y-coil : zero = 2085 
Aaeomlier | 2eCGO ——) 2559 


The voltage values normalized to +5 volts are placed in new 
Mmeamedbbays. At this point the next block of data is read off 
the digital tape and goes through the same process. VOLTR and 
VOLTS analyze eight blocks or 34 minutes of data. The linear 
arrays holding the voltage data are dimensioned to 65,536 
elements given a sampling rate of 32 samples/second. For every 
one of the 65,536 voltage values there is a time value in 
another array. Thus four arrays are required to handle the 
data, three arrays containing voltage data and one that 
contains the respective times for each of the measurements. At 
this point, the VOLTR program calls the NONIMSL subroutine 
DRAWP to plot the data. The VOLTS program uses a smoothing 
algorithm to reduce the amount of high frequency hash. 

ive sneering algorithm is aydouble wunnang sooint 


average executed on all three axes of data. Since I was 
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interested in periods of approximately 10-45 seconds, I 
averaged over 144 points. With a running average over 23 
percent of the oscillation, the unwanted background can be 
removed without destroying the oscillations. Averages greater 
than 50 percent of the period would partially remove the 
oscillation of interest. The construction of the algorithm 
makes changing the number of points averaged over trivial. 

One need only change a loop index and three divisors to obtain 
a different average. The smoothing algorithm is provided as 
maple 3.1. The computer code for both VOLTR and VOLTS is 


attached in Appendix A and Appendix B. 


2. Time Series Voltage Data Analvsis 


One can gain a great deal of insight into the contents 
of the data by looking at the smoothed and unsmoothed voltage 
Beoducts. Figures 3.1-3.3 are unsmoothed voltage plots 
representing the first 34 minutes of the digital data tape GMDT 
ieerecorded 17 August 82, 1301-1408 local time. A one minute 
tape advance was used in producing this data. Immediately one 
notices the large amount of 60 hertz aliasing in the time series 
unsmoothed voltage. The 60 hertz aliasing appears as 4 hertz 
because of the 32 sample/second sampling rate. The sensing 
system has a low pass filter that should be improved to remove 
the aliasing problem. The 4 hertz background can be seen more 
clearly in Figures 3.4-3.6. Above the 60 hertz aliasing 
background, we see voltage spikes ranging in magnitude from l 


volt to over 4 volts. Expanding the time scale to 50 
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Smoothing Algorithm 
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seconds of data, the structure of the voltage spikes becomes 
apparent in Figures 3.4-3.6. The "square well" shape is 
characteristic of PCM dropouts or instances during the PCM to 
digital conversion synchronization is lost between the tape 
recorder containing the PCM tape and the decoding unit. We do 
not know whether PCM dropouts always occur on all three 
channels, but current evidence suggests dropouts are 
Simultanoues on all channels. 

Figures 3./7-3.9 again show the first 34 minutes of GMDT 
ll only smoothed by the 144 point running average. The smooth- 
ing algorithm removes the 60 hertz alising. The PCM dropouts 
appear as "glitches" in the otherwise quasi sinusoid charac- 
teristic of a micropulsation. 

We should note that the unusually large amount of hash 
present between 2 hertz and 10 hertz in the PSD's and the 
average slope of the PSD curve for those frequencies may in part 
be due to the PCM dropouts that occur on the tapes. A 
“deglitching" algorithm could be developed to remove the sudden 
large amplitude structures injected by the system into the 


voltage data if hardware modifications fail. 


B. TIME SERIES MAGNETIC FIELD DATA 
1. Magnetic Field Software 
The time series magnetic field software is basically 
the same as that which produces the PSD's, etc. [Ref. 1 and 
Ref. 3]. The differences begin after the forward Fourier 


transform and the application of the transfer function. We 
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should note that the transfer function as designed ignores 
phase information since it is just a least squares fit-to data 
Tike that of Figure 2.2. Once the transfer function is 
applied, the reverse transform is performed. We obtain the 
complex absolute values of the magnetic field data (in 
nanoteslas) and place the data in new real arrays. To insure 
that there are no discontinuities between blocks of data, 
Figures 3.10a, 3.lla and 3.1l2a, all the blocks are connected to 
the first by adding or subtracting a constant to each element, 
Figures 3.10b, 3.llb, 3.12b. The data now undergoes the same 
144 point double running average as the voltage to bring out 
fluctuations with periods of 10-45 seconds. Again DRAWP is used 
to plot the data. The computer program LVFTC1 produces the 
smoothed magnetic field plots. It is provided in 
Appendix C. 
2. Smoothed Magnetic Field Data Analysis 

Again looking at the first 34 minutes of digital data 
tape GMDT 11, we see in the magnetic field data the same 
exceptional structures as seen in the voltage data, Figures 
3.10-3.12. The glitches, however, do highlight one problem 
rather well. From Figures 3.7-3.9 we note that positive 
voltage spikes produce negative magnetic field spikes. This 
would indicate a pi phase shift is introduced either in the 
Seaiisfer —funcELONS Or in the digital Fourler transform. After 
investigating another digital tape GMDT 3A, 18 August 82, 


0121-0251 local, we see another behavior. Figures 3.13-3.15 
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are the smoothed voltages for the first 34 minutes of data 
on GMDT3A after a one minute tape advance. Figures 3.16-3.18 
are the smoothed magnetic fields for the same time period. The 
X-coil voltages show three positive spikes while the X-coil 
magnetic field shows three negative spikes. This supports the 
theory of a pi phase shift. However, the Y-coil voltage shows 
two positive spikes and one negative spike, while the Y-coil 
magnetic field has three positive spikes. Thus it appears we 
have a rather arbitrary phase problem that may find its origin 
in the transfer functions that do not maintain the correct 
phase. 

From Figures 3.16-3.18 we note exceptional features 
that do not correspond to voltage spikes. The exceptional 
features are separated by time periods of 256 seconds which is 
the size of the blocks analyzed. These features persist with 
the data block connection scheme and 144 point smoothing. Their 
origin is not understood. Some data blocks look like Figure 
3.19 with the end points having about the same value. However, 
some have end points as illustrated in Figure 3.20. The varying 
end point values are also shown in Figures 3.10a, 3.lla and 
3.l12a. This may be a source of the large aoe of hash between 
2 hertz and 10 hertz seen in the PSD's. 

After looking at the magnetic field data it was obvious 
a validation experiment was needed for the sensing system and 


software. 
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migube 3.19. One Frame of Data with End Points Having 
the Same Value 


Figure 3.20. One Frame of Data with End Points Having 
Different Values 





IV. SYSTEM AND SOFTWARE VALIDATION : 


A. EXPERIMENTAL APPARATUS 

The inconsistencies in phase between the time series 
voltage and magnetic field plots prompted an experiment to 
validate the sensing system and software. The sensing 
apparatus consisted of the X-coil and a Schoenstadt fluxgate 
magnetometer. The sensing axis of both sensors were horizontal 
and oriented North-South. Figure 4.1 shows the test set-up. 
Instead of an on axis arrangement, the set-up shown was used 
because of limited resources. The fluxgate magnetometer is 
closer to the source since it has less sensitivity at the 
frequency of interest than the coils. Recorded on the usual 
X-channel was the output from the X-coil. The Y-channel 
contained data from the coil after passing through an additional 
Krohn-Hite model 3321 lowpass filter with a1Ohertz cutoff. On 
the Z-channel was the output from the fluxgate magnetometer. 
All three channels were recorded simultaneously and a 64 
sample/second sample rate was used. An on-site chart recorder 
monitored the fluxgate and double filtered coil voltage output 
real time. The experiment produced two digital data tapes each 
approximately 46 minutes in length denoted GMTT1A and GMTTILB. 

Artificial magnetic fields were introduced into the region 
by applying sinusoidal currents from a Wavetek signal generator 
to a cylinder (diameter = 0.13m) wrapped with 150 turns of 


wire. Location of the source is noted in Figure 4.1. The 
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applied sinusoidal currents varied in frequency from 0.005 

hertz to 10 hertz. For each frequency, two magnitude currents 
were applied to the source to demonstrate a linear response with 
amplitude. 

The fluxgate magnetometer output is correct in phase and 
amplitude to at least 1 hertz. The transfer function between 
voltage and magnetic field for the fluxgate magnetometer is a 
constant of 10 nanoteslas/volt. Thus the fluxgate provides an 
excellent check for the coil transfer functions, the RD data 


stripping subroutine and the digital fourier transform. 


Pee xPERIMENT RESULTS 
1. Voltage 

From the real time chart record of voltage for the 
fluxgate and double filtered coil, we can check the computer 
generated voltages using VOLTR and VOLTS. Table 4.1 shows the 
chart record and computer generated voltages for the fluxgate 
and double filtered coil at various frequencies. Only the 
large amplitude oscillations are measured for the double 
filtered coil because of poor signal-to-noise, for the small 
amplitude oscillations. 

Table 4.1 shows good agreement between the chart record 
and computer output for the double filtered coil. However, 
the fluxgate magnetometer voltages undergo a constant 
Beaicnerptacton Of 0.631 from chart record to computer plot. 


This may be explained by a problem in the data stripping 
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algorithm. The present agreement between the chart record and 
computer generated voltages is shown in Table 4.2. : 

Another more serious problem became obvious during the 
examination of the computer generated voltage plots. For 
approximately the first 20 minutes of data recording the 
fluxgate magnetometer was not connected to the PCM board. 
Figures 4.2-4.4 show the two respective sensing devices. 
Figures 4.2 and 4.3 both contain "real" signal while Figure 
4.4 does not contain "real" signal until the amplitude 
discontinuity indicating connection of the fluxgate 
magnetometer to the PCM board. Figures 4.2-4.4 indicate that 
"pick-up" 1s occurring probably between the channels at the PCM 
board. The amplitude of the cross-talk voltage is about 0.15 
volts peak-to-peak at the frequency of 0.1 hertz. The 
amplitude of the "real" signal on the double filtered coil is 
Mee? volts at O.1 hertz. Thus 51 percent of the “real” signal 
1s coupled to the disconnected channel. Coupling between 
channels of this magnitude raises serious questions about the 
computer generated data, especially for coherence. Discovering 
the amount of the channel pick-up if the channel is loaded 
requires further experimentation. 

Figures 4.5-4.7 show an interesting feature of the coil 
data versus fluxgate magnetometer data. On the channels 
containing the coil output we see noise on top of the signal. 
This noise is not present on the channel containing the 


fluxgate magnetometer data. Since all three channels passed 
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Figure 4.4. 





through the PCM encoder and decoder, the noise is probably in 
the coil amplifiers or the coils themselves. An experiment to 
investigate the true origin of the noise should be done. 

2. Magnetic Field 

Given the dipole equation [Ref. 6] and the experimental 
arrangement in Figure 4.1, the ratio of magnitude of the 
magnetic field seen by the fluxgate to that seen by the coil 
can be calculated. From the geometry in Figure 4.1 the field 
at the coil would be a factor of 8 less than that seen by the 
fluxgate. This is because the source is a dipole and the field 
of a dipole falls off as 1/R . Detailed consideration of the 
geometry gives a ratio of approximately 8.09. Table 4.3 
contains the magnetic field measured at the fluxgate 
magnetometer. Using the ratio calculated from the geometry, the 
expected magnetic field amplitude at the coil is listed in 
Table 4.3. 

A modified form of computer code LFVTC1l generated the 
magnetic field versus time plots for the test tapes. Table 4.3 
shows the computer generated peak-to-peak magnetic field 
values read from the computer plots. From Table 4.3 it appears 
that the coil magnetic field amplitude is in error by a factor 
of approximately three (with exception to the low frequencies). 
One can conclude that a problem exists in the software. 

Figures 4.8-4.9 show the magnetic field data for the 
second seventeen minutes of digital data tape GMTTIA. Figure 


4.9 is the magnetic field measured by the fluxgate magnetometer 
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and Figure 4.8 is the magnetic field seen by the coil. For the 
coil we see discontinuities every 256 seconds corresponding to 
the block size of data. This indicates that the problems in the 
data block end points, as mentioned in Section III, do exist. 

It 1s not clear how to deal with this. 

The first 300 seconds of Figures 4.8-4.9 is expanded in 
Figures 4.11-4.12. With the expanded data a phase difference 
between coil and fluxgate can be determined. By checking the 
phase differece, we can determine whether the phase independent 
meemcrer Lunction is introducing error. For the 0.1 hertz field 
plotted in Figures 4.10-4.11 the coil lags the fluxgate by 
approximately 72 degrees. In Figures 4.12-4.13 we see that at 
O.O0l hertz the coil sensitivity has fallen and a reliable phase 
difference between the fluxgate and coil cannot be made. 

Figure 4.14 shows the relationship between phase angle and 
frequency. 

The above experiment shows the advantages of having a 
meabetime ‘ground truth record. “With this, one can avoid 
“working in the blind". Computer products will have a check. 
Thus, a multi-channel chart recorder should be purchased to 
graph in real time the x, y, and z inputs to the PCM board. In 
addition, electronic circuitry should be developed and 
incorporated so that the chart records represent at least a low 
pass filtering of the coil voltages. Ideally, conversion of 
voltage to magnetic field could occur in real time at this point 


given the appropriate electronics. 
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V. MICROPULSATIONS 


A. THEORY 

Currently there is no theory that completely models the 
observed micropulsations and geomagnetic noise. It is 
generally accepted that the generating mechanism for 
micropulsations is an interaction between the solar wind and 
the earth's magnetosphere [Ref. 5]. The observed 
micropulsations are classified by their period and regularity 
LRef. 5]. Table 5.1 shows the micropulsation classification. 
The smoothing algorithm applied to the magnetic field data is 
designed to highlight micropulsations with periods between 10-45 
seconds. Specifically we are interested in micropulsations 
classified as Pe 3. Pe 3 micropulsations are generally regular 
and have a more or less distinct period. Thus in the PSD plots 
we will be able to identify the micropulsation by an apparent 
anomaly at its frequency. Micropulsations have amplitudes of 
approximately 0.5 nT. This small amplitude is very small 
compared to a main magnetic field value of 50,000 nT, hence the 
term "micro" pulsations. Geomagnetic fluctuations are a 
combination of micropulsations and a general background noise. 
Micropulsation events grow out of the ever present background. 
The time series magnetic field plots will provide a means of 


obtaining the frequency of the micropulsation. 


Aro 





Hable wo.) Mieropulscation Classification 


Notation Period (seconds) 
Bol O25 
Pez 5-10 
Po3 10-45 
Pc4 45-150 
Pc5 150-600 


alsaleet 





B. MICROPULSATION DATA ANALYSIS 

The first 1200 seconds of data after a 60 second tape 
advance on digital data tape GMDT 11, 17 August 82, 1301-1408 
Local, may contain a micropulsation event. Figures 3.10-3.12 
show the first 34 minutes of GMDT 11 after a 60 second advance. 
To determine whether the regular pulsations we see in the first 
third of Figures 3.10-3.12 (disregarding the suspicious "cusps") 
is a micropulsation, the scale must be expanded. Figures 5.1- 
5.3 show the regular pulsations on an expanded scale. From 
Baoures 5.1-5.3 the period of the oscillation is calculated to 
be 17.2 seconds. Thus this section of data contains a type Pc 3 
micropulsation. 

Executing the computer code LANDXYZ [Ref. 1] on this 
section of data, we generate the PSD, coherence, degree of 
Semarrzation, and ellipticity. Table 5.2 contains the 
amplitude, frequency, coherence, etc. for the micropulsation. 

At the frecuency of the micropulsation one would expect a 
good coherence, degree of polarization and ellipticity since 
the coil system is oriented arbitrarily with each coil sensing 
the micropulsation. The observed micropulsation exhibits a 
high coherence as shown in Table 5.2. The observed micropulsa- 
tion has a high degree of polarization. Micropulsations are 
generally considered to be left-hand elliptically polarized 
[Ref. 5]. The observed micropulsation is elliptically 
polarized. 

If the errors in phase and amplitude that were shown to 


exist in Section IV are constant, then that should not be a 
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Frequency (Hz): 0.058 


Micropulsation Properties 
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Peak-to-Peak Amplitude (nT): Amplitude = ae v7, = Oro 7 
PSD (dB) (Ref lnT /y,): X-Coil: -34dB 
Y=Coil: -4qdB 
Z=COl es SaB 
Coherence: X=-Y Coils: 0.99 
Y=Z Geil ss 0.99 
f= he On See Oo 
Wegree of Polarization: X-Y Plane: 90.99 
Y-Z Plane: 0.99 
Z-K- Plane: 0299 
Epeieeot city: X-Y Plane: 0.45 (left hand) 
Y-Z Plane: 0.10 (right hand) 
Z-X Plane: 0.30 (left hand) 
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serious problem since we are dealing with a single frequency. 
Thus, the observed micropulsation exhibited definite structure 


decidedly different than the background. 


ieee 





VI. CONCLUSIONS 


A type PC-3 micropulsation event of peak-to-peak amplitude 
(0.37 nT) and period (17.2 seconds) was identified in a section 
of data. At the frequency of the micropulsation, the coherence 
between the three orthogonal coil sensors was about 0.99 for 
miemxy~-pair, 0.99 for the xzZ-pair, and 0.99 for the yzZ-pair. 
The degree of polarization was 0.99. The ellipticity was about 
Oe 0. 

The origin of the regular "hump"-like structure in the 
PSD's of previous coil data was found to be artificial and not 
representative in any way of natural phenomena in the 
geomagnetic field. Previous conclusions based on a frequency 
domain analysis of coil data characterized by a regular "“hump' 
structure (strong or weak) in the PSD plots are known to be 
unsubstantiated. 

Compared to a PSD frequency domain analysis of the data 
set, an analysis of geomagnetic field fluctuations in the time 
domain easily identifies data sets containing a micropulsation 
event. A double running average routine developed and applied 
to the time series data sets was found to be helpful in 
nes itonting the micropulsation. 

Detailed conclusions about the geomagnetic field from the 
coil measurements require a balance of analysis in the 


frequency and time domain. 


is 





Pumedaeccurate reproduction Of a time series plot for the 
geomagnetic field oscillations will require modification to the 
existing software. 

Comparisons in the time or frequency domain between any 
pair of coils in the orthogonal system might be misleading as a 
result of cross talk between electrical inputs. The suspicion 
of cross talk is based on the observed vick-up on an unloaded 
input of the contents of an adjacent loaded and active input. 

Further micropulsation studies should continue. This 
initial investigation of the coherence, polarization, and 
ellipticity features of micropulsation events shows promise 
that with more data, one might suggest novel physical 


mechanisms and models for the origin of micropulsations. 


To 





APPENDIX A 


VOLTR COMPUTER PROGRAM 


Appendix A contains the VOLTR computer program that 


provides unsmoothed time series voltage plots. VOLTR uses the 


subroutine RD to strip voltage data off digital tapes. The 


integer voltage data is normalized to + 5V and plotted using 


the DRAWP subroutine. 
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APPENDIX B 


VOLTS COMPUTER PROGRAM 


Appendix B contains the VOLTS computer program that 
provides smoothed time series voltage plots. VOLTS uses the 
Subroutine RD to strip voltage data off digital tapes. The 
integer voltage data is normalized to + 5V. A double 144 point 
running average is applied to the data. The DRAWP subroutine 


plots the data. 
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PEER WDEK CC 


LFVTC1 COMPUTER PROGRAM 


Appendix C contains the LFVTCi computer program that 
provides smoothed time series magnetic field plots. LFVTCl 
uses the subroutine RD to strip voltage data off digital tapes. 
Maesvoltage data 15 normalized to + 5V. The voltage data is 
fast fourier transformed using the FOURT subroutine. The 
transfer function is applied and forward transform taken. The 
resulting magnetic field data undergoes a double 144 point 


running average and is plotted using the DRAWP subroutine. 
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